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ABSTRACT 



The reduction of integral-field spectrograph (IFS) data is demanding work. Many repetitive operations are required in order to convert 
raw data into, typically a large number of, spectra. This effort can be markedly simplified through the use of a tool or pipeline, which is 
designed to complete many of the repetitive operations without human interaction. Here we present our semi-automatic data-reduction 
tool p3d that is designed to be used with fiber-fed IFSs. Important components of p3d include a novel algorithm for automatic finding 
and tracing of spectra on the detector, and two methods of optimal spectrum extraction in addition to standard aperture extraction. p3d 
also provides tools to combine several images, perform wavelength calibration and fiat field data. p3d is at the moment configured 
for four IFSs. In order to evaluate its performance we have tested the different components of the tool. For these tests we used both 
simulated and observational data. We demonstrate that for three of the IFSs a correction for so-called cross-talk due to overlapping 
spectra on the detector is required. Without such a correction spectra will be inaccurate, in particular if there is a significant intensity 
gradient across the object. Our tests showed that p3d is able to produce accurate results. p3d is a highly general and freely available 
tool. It is easily extended to include improved algorithms, new visualization tools and support for additional instruments. The program 
code can be downloaded from the p3D-project web site http://p3d.sourceforge.net. 
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1. Introduction 

With integral integral field spectrographs (IFSs) an extended 
area on the sky can be spectroscopically mapped, under the same 
observing conditions, in one single exposure. In order to fit all 
simultaneously observed spectra onto the detector, the field-of- 
view of the integral field unit (IFU) that provides the spatial sam- 
pling on the sky is relatively small. The footprint of IFUs typi- 
cally ranges from several arcsec, e.g. PMAS, to about one ar- 
cmin for VIMOS. In an alternative configuration, IFUs are used 
to map a much larger area, but with sparse sampling, e.g. PPAK 
and VIRUS-P. In Table Q] we list several of the existing fiber- 
fed IFUs, the telescope they are mounted on, an d the names of 
the corresponding pipelines. See lBershadvl d2009t) for a complete 
list. Regardless of the actual IFS configuration, the raw data of 
fiber-fed IFSs always consist of hundreds, or even thousands, of 
spectra per exposure. The reduction of these data consists of pro- 
cessing each spectrum individually and is therefore highly repet- 
itive work. A general data-reduction tool is highly desirable, that 
automates the reduction steps, yet allows the user to interactively 
inspect and optimize parameters when required. The purpose of 
p3d is to provide such capabilities and thereby facilitate the sci- 
entific exploitation of IFSs. 

Special purpose data-reduction pipelines exist for most IFUs. 
Two additional reduction packages that are more general in 
their functionality than the ones listed in Table [1] and are suit- 



able to use with fiber-fed IFUs, are r3d (ISanchezH2006l here- 
after S06) and iraiQ. These tools, however, require significant 
amounts of time-consuming manual interaction in order to get 
the best out of them. While developed initially for the PMAS 
spectrograph, the data-reduction tasks of p3d (and p3d_online) 
are very generally formulated and are equally applicable to data 
obtained with other fiber-fed IFSs. The suitability of p3d for sev- 
eral other IFSs is demonstrated with a variety of scientific re- 
sults at an early stage, when no such tools were available yet 
for those instruments shortly after commissioning, e.g. MP FS 
dRoth et al.ll2004t iLehmann et al J 120051: iFabrika et alj|2005h or 



VIMOS (iMonreal-Ibero et alJl2005HVillar-Martm et al.ll2006l) . 

In this paper we present a generalized v ersion of p3d, th at is a 
complete rewrite of the previous version of lBeck er (2002). This 
general IFS reduction tool includes support for PMAS/LARR, 
PMAS/PPAK, VIRUS-P and SPIRAL, and can be readily ex- 
tended to additional IFUs. Key features of p3d are: 

- A single, freely available, and easy to install, program pack- 
age with support for several IFSs and computing platforms. 

- Calculation and propagation of errors through all steps. 

- The option to choose between aperture extraction and two 
methods of optimal extraction. 

- A possibility to store information about all performed oper- 
ations in a log file. 

- Interactive and integrated visualization tools to allow the 
user to examine intermediate and final products. 



Send offprint requests to: CSandin@aip.de 

* Based in part on observations collected at the Centra Astronomico 
Hispano Aleman (CAHA), operated jointly by the Max-Planck Institut 
fiir Astronomie and the Instituto de Astroffsica de Andalucia (CSIC). 



1 The Image Reduction and Analysis Facility iraf is distributed by 
the National Optical Astronomy Observatories which is operated by the 
association of Universities for Research in Astronomy, Inc. under coop- 
erative agreement with the National Science Foundation. 
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Table 1. A list of fiber-fed IFUs and their respective data-reduction pipelines 
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PPAK 
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SPIRAL 
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VLT/UT3 
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VIMOS-IFU 
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VIPGI 

VTMOS PIPELINE 


8a 


McDonald 2.7m 


VIRUS-P 


VIRUS-P 
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Notes. The name of the telescope, the spectrograph and the IFU are given in Cols. 1-3. In Col. 4 we specify the number of detectors of the IFU. 
Column 5 specifies the main instrument reference paper, and Cols. 6-7 give the name and reference of instrument-specific reduction tools/pipelines. 

References. lAvila et al l d2003h. '{B~lecha et all d2000h. '*|Pa"sauini et alj (120001): ^Allington-Smith et alJ d2002h: ISchmo 1 et all d2004. 
3 lBolton & Buries! d2007h: jArribas et al.l dl998T> : iRoth et all d2005[). 5i lBeckeri d2002l):"lK'elz et al.l J2006h : ISmith et al.l d2004l) . 7 lSharp et al.l 
d2006h : iLeFevre et all d2003h . IScodeggio et all d2005h : Hffill et alJ d2008h . lAdams et all d201Ch . 



The code, furthermore, includes full run-time error handling and 
both the code and supplementary data files are fully commented. 

This paper is laid out as follows. In Sect. [2] we first describe 
the goals and setup of p3d. The data-reduction algorithms and 
their implementation are thereafter detailed in Sect. [3] We dis- 
cuss and analyze the outcome of the tool, and make a compari- 
son with corresponding outcome of iraf, for some of the tasks, 
in Sect. [4] Finally, we close the paper with our conclusions in 
Sect.|5] 

2. About the objective and setup of p3d 

Our goal with p3d is to provide a general, flexible, fast and reli- 
able tool for reduction of fiber-fed IFU data. p3d is not intended 
to be a tool that handles every possible task from data reduction 
to data analysis. Instead it focuses on the tedious and repetitive 
tasks, which are required to convert raw data into wavelength- 
calibrated spectra, cf. Sect. [3] We also want to provide a user- 
friendly tool where the required input from the user is kept to a 
minimum. It is recommended, however, that the user has a basic 
knowledge of the data-reduction process, and the related numer- 
ical problems and methods, in order to allow a full exploitation 
of all benefits of p3d. 

T he collection of routines that m ake up p3d were from the 
start dBeckerfl2002l iRoth et al.ll2005l) written in the Interactive 
Data Languag^3 (IDL). Supplementary routines, which are 
mainly related to file-IO, are used from the public ly available 
toolkit astro-lib of nasaB The mpfit routine of iMarkwardtl 
(l2009lFl is also used when fitting line profiles (for the optimal 
spectrum extraction). The current version of p3d - that is a public 
release under GPL-v3 - is a complete rewrite of earlier versions 
and remains a graphical tool (GUI), although reduction tasks can 
also be completed without the GUI. We re-designed the tool to 
work with all platforms which are supported by IDL0 (version 
6.2 or higher). We stress that p3d can be used with full function- 
ality without an IDL license; using the IDL Virtual Machine, 
all required scripts are provided to use this functionality (with 

2 http://www.ittvis.com 

3 http://idlastro.gsfc.nasa.gov 

4 http://purl.com/net/mpfit 

5 p3d uses widgets abundantly and can therefore not be used with 
GDL (cf. http://gnudatalanguage.sourceforge.net). 



all UNIX-type platforms). The installation procedure consists of 
setting up IDL to include the files of p3d, and the supplemen- 
tary routines, in the path of IDL. The software, updates, docu- 
mentation and tutorials can be found at the p3D-project web site, 
http://p3d.sourceforge.net. 

All routines and their input and output parameters are doc- 
umented. When p3d is run using the GUI the functionality of 
every interactive part is described with a comment in a sta- 
tus indication field, and information on all activity is option- 
ally saved to a log file. We minimized the run-time by replac- 
ing time-consuming FOR-loops with intrinsic IDL functions. All 
parts of the GUI can be used with screen sizes from 1024x600 
and larger. Moreover, the program, currently, consists of about 
45 000 lines of code. Any redundancy is kept to a minimum 
by adhering to the "Don't Repeat Yourself" principle of cod- 
ing dHunt & Thomasl 19991) . Our main priority is to correct algo- 
rithm issues, which lead to erroneous scientific output, as soon as 
possible. Avoidable issues, which are usually related to different 
parts of the GUI, are corrected as our time permits. 

Currently p3d is configured for four IFUs: the lens array 
(LARR) and PPAK IFUs of the PMAS instrumen t (for both 
the o l d 2k x 4k CCD and the new 4k x 4k CCD, iKelz et al.l 
2006} IRoth et al.ll2005h that is mounted on the 3.5 m telescope 
at Calar Alto, t he VIRUS-P IFU (for both bundle 1 and the 
newer bundle 2. iHill et al.ll2008l) that is mounted on the 2.7 m 
Harlan J. Smith telesc ope at the McDona ld Observatory, and the 
SPIRAL IFU (see e.g. ISmith et al]|2004l) that is mounted on the 
3.9 m Anglo-Australian telescope. Support for additional instru- 
ments can be added through appropriately formatted configura- 
tion files. 

3. Components of the data reduction 

p3d contains a set of routines to do the following five tasks: cre- 
ate a master bias, trace all spectra on the detector, create a disper- 
sion mask, create a flat field, and extract spectra in object data. 
We now describe these tasks in more detail. 

For each task several raw data images can be combined by 
p3d in order to increase the signal-to-noise and remove cosmic 
rays. The default combination method is to use a stack of at least 
three images, where the minimum and maximum values of every 
pixel are thrown away before calculating a (min/max-)average 
of all images. Optionally an average or a median can be used. 
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Fig. 1. This figure shows two adjacent PPAK raw data sec- 
tions highlighting one group of 32 spectra. We show continuum 
lamp data on the left-hand side and object data, about Ha, on 
the right-hand side. Darker shades indicate stronger intensities. 
Horizontal tick marks, which we show on the left hand side of 
the image, indicate positions of individual spectra. The vertical 
dashed line indicates the search direction of step 1 of the tracing 
algorithm. The horizontal dashed line indicates, for one spec- 
trum, the trace that results after step 2. For further details, see 
Sect. El 

Prescan- and overscan-regions in the raw data are only removed 
immediately before spectra are extracted; this must be accounted 
for if external tools are used to create intermediate products. 
Spectrum extraction can be done using either aperture extrac- 
tion, our own so-called modified optimal extraction, or using 
a multi-profile deconvolution optimal extraction (cf. Sect. 13.21 ). 
Both optimal extraction methods are able to correct for so-called 
cross-talk due to overlapping spectra. Currently p3d does not in- 
clude tasks to apply a bad pixel map, remove scattered light, 
subtract the skjQ, perform flux calibration, compose dithered or 
mosaiced combined frames, or any kind of binning of spectra. 

In the first step a master bias image is created by combining 
a set of at least three bias images. This master bias is subtracted 
from the raw data in all consecutive steps. 

In the second step the position of every spectrum is deter- 
mined along the dispersion axis with a well-illuminated calibra- 
tion exposure, using a continuum lamp or twilight flats. The re- 
sulting trace mask is used in all consecutive steps when spectra 
are extracted. Finally, cross-dispersion profiles are calculated for 
all spectra and wavelengths in order to allow optimal extraction. 

In the third, optional, step a dispersion correction is deter- 
mined for every spectrum, using one or several arc lamp expo- 
sures. An extracted flat field image is created in the fourth, also 
optional, step in order to correct for wavelength-dependent vari- 
ations of every spectrum and for differences in the fiber-to-fiber 
throughput. In the fifth, and final, step all spectra are extracted 
from object exposures, optionally applying first the dispersion 
correction and thereafter the flat-field correction. 

In order to illustrate typical properties of IFU raw data we 
show a section of object data in Fig. [1] together with the cor- 
responding continuum lamp data, which were taken with the 
PPAK IFU. Every horizontal line marks the position of a spec- 
trum, the two brightest lines (that are only seen in the left hand 
side image) show calibration fiber spectra of PPAK (these are 
not used in p3d). The strong variation of emission line intensi- 
ties of different spectra (on the cross-dispersion axis ) is due to 
the spatially irregular mapping of fibers in PPAK (cf. Kelz et al] 



120061) . The two data sets that are shown next to each other il- 
lustrate their interdependence - spectrum positions in the object 
data should match those in the continuum (calibration) data well 
in order to extract spectra properly. 

In the following we assign all variable parameters a symbol 
so that they can be found easily inside the program code. For 
a quick reference we collected all symbols, and their respective 
default value, in Table [2] the meaning of the parameters are in- 
troduced step-by-step below. The CCD readout binning parame- 
ters are, furthermore, denoted by bin,i (dispersion axis) and bin ) 
(cross-dispersion axis); these parameters are set to 1 or 2. 

Next we describe how p3d handles the separate data- 
reduction tasks. Our novel spectrum tracing algorithm is de- 
scribed in Sect. 13.11 Thereafter we describe our approach to 
spectrum extraction in Sect. 13.21 how we wavelength calibrate 
the data in Sect. 13.31 how we flat field data in Sect. 13.41 and how 
object spectra are extracted in Sect. 13.51 



6 There are several ways to subtract the contribution of sky emission 
lines from object data. For the initial release of p3d we have chosen 
not to include any sky subtraction as it is difficult to define a general 
procedure for it. 



3.1. An automated method for finding and tracing spectra in 
IFU raw data 

Before spectra can be extracted they must be found and their ex- 
tent be traced along the dispersion axis on the detector. In order 
to allow accurate tracing a continuum lamp exposure should be 
used where spectra are clearly visible across the full wavelength 
range. Most IFU-instruments are, moreover, affected by flexure, 
which results in changing positions of spectra on the detector as 
the telescope moves over time; examples of instruments without 
such flexure are VIRUS -P and INTEGRAL. If the flexure is sig- 
nificant it is important to trace with a continuum lamp exposure, 
which is sampled close in time to the object exposure, in order 
to extract data from the correct region on the detector. It may 
also be an option to calculate one trace mask and then shift it 
using knowledge about shifted positions of selected lines of arc 
images. 

We here present an automatic algorithm for spectrum tracing 
that involves several steps to assure that both all expected 
spectra of an instrument are found and traced accurately, and 
that no noise feature is kept as a spectrum. The method is 
general and should require little if no modifications once it is 
setup for an instrument. No assumptions are made regarding 
exact positions of spectra on the detector. Instead information 
is required about their expected number, their orientation, and 
their ordering. Cosmic ray hits could, if numerous, have some 
effect on the outcome; they should in that case be removed in 
advance. 

The algorithm is split into two steps. At first all n spectra are 
searched along the cross-dispersion (spatial) axis for one pixel, 
or a set of averaged pixels, on the dispersion axis. Thereafter they 
are traced along the dispersion axis for the remaining pixels. All 
instrument-specific default parameter values of this procedure 
are shown in Table [2] for all considered IFUs. We next describe 
the two tracing steps separately. 



3.1 .1 . Step 1 : finding the spectrum positions 

At first local maxima, representing positions of individual spec- 
tra, are located on the cross-dispersion axis; in this process a set 
of adjacent pixels is used on the dispersion axis, which is typi- 
cally centered on the middle pixel of the CCD. It is assumed that 
spectra are fairly well aligned with either axis on the detector. 
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Table 2. General and instrument-specific parameters of the data-reduction algorithms of p3d 



Step 


Property All 


PMAS 


VIRUS-P 


SPIRAL 


Parameter name in p3d 






LARR 


PPAK 


bundle 1/bundle 2 


both arms 




Automatic tracing algorithm, cf. Sect. 13. 11 












n 


256 


382 


Z40 


^ 1 1 

J tz 


SPNUM 


la 


\n„ — 1 ] bin 1 11 


40 


40 


29 


30 


FINDWIDTH.TR 






0.9 


0.9 


0.9 


0.95 


CUT_TR 


lb 


d bin t 


12.5 


9.5 


8.0 


4.0 


DIST_TR 




(5 mln bin + 


4.0 


2.5 


2.0 


1.0 


DMIN_TR 




<Wbin t 


5.5 


4.0 


4.0 


2.5 




lc 


(iu cc - l)bin-| /2 2 










CENTERVAR_TR 




G w bin t 


4.0 


4.5 


J.U/4.Z 


Z.J 


FWHM_TR 




Hit 9 










IN 1 1 UK Al _1K 


2 


f 10 
(u>«-l)biiu/2 10 
(ic c - 1) /2 10 
(Wj-D/2 5 










REFINDDIST_TR 
REFINDWIDTH_TR 
SMOWIDTH.TR 
DISPSMOWIDTH_TR 


Spectrum extraction, cf. Sects. 13. 2. 1113. 2. 31 












(jr w - l)bin t /2 


6.0 


4.5 


3.4/3.3 


1.5 


PROFWIDTH_EX 




(x w - l)bin t /2 


8.0 


8.0 


8.0 


5.0 


PROFWIDTH_CTEX 




^mpd 










MPDNPROF 


Line profile fitting, cf. Sects. 13. 2. 2113. 2. 31 












OJ/ 


16 


31-32 


21-44/14-29 


32 


LPROFN_TR 




« f bin,! 


42 


12 


104 


22 


LPROFDINT_TR 



Dispersion mask creation, cf. Sect. [33 
P 4 
W - l)bin, i /2 4 

Flat fielding, cf. Sect.[23 
(«%-l)bin,,/2 20 
Pft 7 



POLYNOMIALORDER.DM 
LINEWIDTH_DM 



SMOWIDTH_FF 
DEG_POLYFIT_FF 



Notes. Excepting the dimensionless properties n, n a , mi, «f, p and p^, the unit is pixels. In the third column we specify default values of general 
parameters. In the last column we give the name as it is specified in the instrument-specific parameter file of the program code. 



The following condition must hold for any local maximum, 

/(/// 1) </(//.) -'/(y-:) • 

where f{y) is the intensity distribution across the cross- 
dispersion axis {y) for every pixel j, and j is a pixel with a local 
maximum. In order to reduce effects of noise, possible cosmic 
ray events, and variations in the detector sensitivity across its 
surface, the local maxima are searched in a set of n p adjacent 
pixels on the dispersion axis (x). For positively identified spec- 
tra we require that the same pixel position in y is found in at least 
a fraction <f of the set of n p pixels in x. By allowing maxima to 
be present also in the next pixel (yj+i) slightly tilted spectra are 
found as well. The statistical probability for spurious detections 
of noise features as spectra, which depends on the instrument 
setup, is nevertheless significant using only this step. 

Next, in step lb, knowledge of the spectrum pattern separa- 
tion is used to filter out spectra from the sequence of local max- 
ima of step la. The main assumption is that consecutive spec- 
tra are separated by an instrument-specific distance (pitch) d(y). 
The first maximum in the sequence is used as a starting point 
of a sequence. In order to match a sequence the distance to the 
subsequent maximum must fulfill 

mxd(y) - 6 min < mxd(y) < mxd(y) + 5 max , 

where mxd is an integer multiple of d (m > 1) and <5 m j n and 
#max specify the permitted deviations. Mismatched maxima be- 
gin another sequence. By allowing such deviations it is possible 
to permit gaps between groups of spectra where the separation is 



non-constant. The spectrum pitch d, (5 m ; n and £ max must be deter- 
mined for every instrument in advance. The spectrum separation 
d(y) is in general constant with y, althou gh with the INTEGR AL 
IFU it varies across the detector (see e.g. lArribas et alJl998l this 
is not yet handled by p3d). In Fig. [2] we show how the separa- 
tion of consecutive spectra varies for the four supported IFUs 
(similar figures can be used to fine-tune d, S m i n and <5 max for ad- 
ditional IFUs). Large separations indicate gaps between groups 
(or banks) of spectra. The separation between spectra belonging 
to separate groups is roughly constant for all four IFUs. Once 
all maxima have been traversed the longest resulting sequence is 
selected as the sequence of real spectra. 

In step lea cross-correlation is made between the sequence 
of maxima of step lb and a pre-defined instrumental pattern, that 
specifies the expected separation and number of spectra. This 
pattern is defined as a list of spectrum gaps, assuming spectra 
are separated as is described in step lb (setting m = 1). For ex- 
ample, if two groups of spectra are separated by a distance cor- 
responding to 2xd then this corresponds to one entry in the list 
of gaps. After the cross-correlation positions of missing spectra, 
due to e.g. dead or unused fibers, are inserted separately. Those 
positions, which could not be identified in steps 1 a and lb, are in- 
terpolated or extrapolated from the position of the nearest found 
spectrum. In this way the number of spectrum positions in the 
returned sequence is always as expected, viz. n (this is a use- 
ful property when working with data of IFUs such as VIMOS 
that has many fibers with poor throughput). After this third step 
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Fig. 2. In this figure we illustrate how adjacent spectra are separated for four IFUs: a) PMAS-LARR (using bin-|- = 2), b) PMAS- 
PPAK (bin t = 2), c) SPIRAL (bin t = 1), and d) VIRUS-P (bundle 1; bin t = 1). The ordinates are logarithmic and common to all four 
panels. The separation to the previous spectrum is indicated with the symbol x. In panel c (d) the lowermost x indicates (the two 
lowermost x indicate) the calculated offset of a dead fiber. For further details, see Sect. 13.1.11 



the probability that a noise feature is identified as a spectrum is 
negligible. 

Finally, more accurate positions of the identified sequence 
of spectra are calculated by weighting with the cross-dispersion 
profile of the data. In this weighting the width of the used 
spectrum section is w cc pixel. The full-width at half-maximum 
(FWHM) of the Gaussian profile that is used in the weighting is 
G w pixel. This procedure is in every case iterated n lt times. 

In Fig. [3] we show an example of the outcome of step 1 of the 
tracing algorithm for calibration data, which were sampled with 
the PPAK IFU. For a direct comparison we show the same sec- 
tion of the spectrum that we present in Fig.[T] Note that 35 out 
of 36 visible spectra were identified in steps la and lb. The po- 
sition of spectrum 225 is hidden in the much stronger spectrum 
224. The position of this spectrum was determined in step lc, 
where the spectrum pattern defined its location (as d pixels dis- 
tant from spectrum 224). All spectra, but those of the calibration 
fibers 192 and 224, are to a lesser degree affected by cross-talk, 
cf. Sect. g2] 

3.1.2. Step 2: tracing along the dispersion axis 

In a second step we calculate spectrum positions for all wave- 
length bins, starting with the resulting weighted positions of step 
1. If spectra are traced individually along the dispersion axis 
there is a risk that some spectra could be lost, and because of 
effects of pixel subsampling an oscillating pattern of individual 
spectra might be found in the trace mask. The amplitude of the 
oscillations decreases with the profile width and becomes signif- 
icant with smaller widths. The associated oscillation period de- 
pends on the angle between the spectrum and the row of pixels, 
and corresponds to the number of pixels in a row of pixels host- 
ing a maximum. If positions are smoothed, using appropriately 
selected smoothing parameters, such oscillations can be made 
smaller, cf. Sect. [471] 

We use the following approach. At first the data set is re- 
sized, by a defined factor /, on the dispersion axis, in order to 
reduce effects of noise. In this process the value at every position 
is averaged over w a pixel. For every bin on the dispersion axis 
we then calculate weighted spectrum positions, starting with the 
already known position of a neighbor wavelength bin. In order 
to remove subpixel-sampling effects the positions are thereafter 
smoothed twice using a box-car of width w c pixel, first using a 
median and then using an average. After all positions are cal- 



culated we smooth the positions on the dispersion axis, using a 
box-car of width u>a pixel. 

Finally we calculate spectrum positions for all pixels on 
the dispersion axis. Using the smoothed and rebinned data 
set remaining positions are calculated by linear interpolation. 
Values at either end are extrapolated. 



3.1 .3. Adjusting the automatic tracing procedure 

The philosophy of the described tracing algorithm is that it is 
automatic, and no modifications of the default values of Table [2] 
should be required. However, in some circumstances it may turn 
out that the algorithm is unable to correctly identify all spectra. 
If this were the case the recommended procedure is to first vary 
£ to see if the number of spectrum misfits can be made smaller 
(step la). A second option is to vary 5 mdx , and maybe also 6^, 
in order to find outliers (step lb). If the number of misfits is 
large it may be necessary to add or remove one entry in the list 
of spectrum gaps (step lc). If the spectrum separation parameter 
d is changed it is in any case necessary to modify this list. 

3.2. Introducing the three spectrum extraction methods 

The flux in every wavelength bin, of any spectrum on the detec- 
tor, is distributed in a profile, or aperture, on the cross-dispersion 
axis. The shape and extent of the profile depends on the instru- 
mental setup. If many spectra are squeezed onto the surface of 
the detector, as is often the case with IFSs, there is likely some 
overlap between profiles of adjacent spectra, resulting in cross- 
talk. The smaller the spectrum separation (pitch) and the larger 
the spectrum width are, the greater the overlap is. By selecting 
an accurate method of spectrum extraction it is mostly possible 
to separate overlapping spectra well and attain both accurate and 
precise values of the flux. Such an approach can, however, be 
very time-consuming, which is why it is worth testing if a sim- 
pler method suffices. 

From here on we assume that the data set o, that contains 
spectra, is bias-subtracted . The varianc e of the flux for every 
pixel k in o is then (cf. e.g. [Howell 20Q(J, 

V , k = \d k -b k \/g k + rl + V h , k , V b , k = r 2 k /n h , (1) 

where d is the raw data with spectra, b is the master bias, g is the 
CCD gain (e~/ADU; ADU is the analog-to-digital unit), r is the 
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Fig. 3. This figure shows the result of step 1 of the tracing algorithm for continuum lamp data of the PPAK IFU, cf. Sect. 13.1.11 In 
this case the used CCD binning setup was bin; xbin-|- = 2x2. The intensity is shown for the centermost pixel on the dispersion axis. 
Only 36 of 382 spectrum positions of PPAK are shown, every identified position is indicated with a vertical line and a spectrum 
number. Spectrum positions found in step la are indicated with a plus symbol (+). Positions found in step lb are indicated with a 
times symbol (x), and finally adopted spectrum positions are indicated with asterisk symbols (*). Note that the two spectra 192 and 
224 are spectra of calibration fibers. 



readout noise (ADU), Vb is the variance of the averaged master 
bias, and n\, is the number of individual bias images, which were 
used when combining the images. A pixel index is added to the 
gain and the readout noise in order to account for instruments 
where these properties vary across the CCD surface (such as is 
the case with the new 4kx4k CCD of PMAS that is read out in 
four blocks, which have to be combined before the data is used). 

Next we describe the standard aperture extraction in 
Sect. 13.2.11 Thereafter we introduce our modified optimal ex- 
traction (MOX) method in Sect. 13.2.21 and how we implement 
the multi -profile deconvolution (MPD) method in Sect. 13.2.31 

3.2.1. Standard aperture extraction 

In our first approach we use an aperture of pre-defined width 
with a top-hat (box-car) profile that is set to 1 within the aper- 
ture, and outside. This method is also ref erred to as tramline 
extraction (see e.g. ISharp & Birchalll FioTol hereafter S10). We 
account for aperture boundaries inside of pixels by adding frac- 
tions of flux of such pixels. The integrated flux / and the cor- 
responding variance V s of every spectrum i and the contributing 
pixels of every profile, for the wavelength bin A, are, 

fu = 2 d,M - b#, V 8 ,tf = 2] v o,M- (2) 

M AO 

We present the instrument-specific default aperture widths, x w , 
in Table[2]for all four IFUs. 

Although almost all flux can be collected with an aperture 
extraction for IFUs such as LARR this is still an inefficient 
method if fluxes are low. If this is the case the contribution of 
readout-noise from the outer pixels of an aperture can become 
significant, and even dominate flux errors, cf. Sect. 14.41 Since 
this is the fastest method of spectrum extraction it is, neverthe- 
less, the default method of p3d. 



3.2.2. Modified optimal extraction (MOX) 

Optimal extraction ( lHornd[T986l hereafter H86) is a more ac- 
curate method than aperture extraction. In this method line pro- 
files are used to weight flux from separate pixels across aper- 
tures. With such line profiles it is also easier to filter out pixels, 
which are hit by cosmic rays. Although the process of extraction 
is straightforward the calculation of line profiles is demanding 
with IFUs. The most direct approach requires a separate profile 
to be calculated for every spectrum and wavelength bin. This is 
not only a time-consuming process, but also makes it difficult 
to correct for effects of noise and cross-talk; S06 attempts such 
an approach and draws a similar conclusion. The number of free 
parameters may be almost equal to the number of pixels, and 
possible effects of pixel subsampling (see above) are neglected. 

Starting at one wavelength we simultaneously fit a group of 
line profiles, using a pre-selected function. Doing this we use the 
same continuum image which we used to calculate traces. The 
number of spectra mi in every group I depends on the instrument 
setup. We show the ranges of m; which we use with the four 
IFUs in Table |2 The function can be selected to be a Gaussian 
function, a Lorentzian function, an approximative Voigt profile, 
or a double Gaussian function. As we currently ignore any broad 
component due to scattered light it appears sufficient to use a 
single Gaussian profile for all IFUs. 

Since the profile width typically changes slightly across the 
detector it seems inappropriate to fit all spectra simultaneously, 
and such an approach is also computationally much more time- 
consuming. For every group of spectra we (by default) assume 
fixed profile center positions using the trace mask (the center po- 
sitions can be fitted as well, if necessary). Using a Gaussian line 
profile the result of the fit consists of the spectrum width G Wi ;, / 
intensities, and the zero-level and the gradient of the linear fit to 
the background. Thereafter the same procedure is carried out for 
a set of additional wavelength bins, and as a last step we interpo- 
late the profile parameters for all intermediate wavelength bins 
using cubic splines. The number of wavelength bins, that are 
used to calculate the fits, tif, should be selected to allow a rea- 



C. Sandin et al.: p3d: a general data-reduction tool for fiber-fed IFUs 



7 



sonable interpolation for intermediate wavelengths using such 
splines. Finally, before spectrum extraction the profile parame- 
ters are used to calculate pixel-based line profiles for every spec- 
trum and wavelength bin. Every profile is thereafter normalized. 

In comparison to "plain" optimal extraction our modified ap- 
proach allows a correction for overlapping spectra (i.e. cross- 
talk). Including masking of pixels, which are hit by cosmic rays, 
and correction for cross-talk, the equation of the modified opti- 
mally extracted flux becomes, 



fu 



£j (i) M ijA p ijA (a jA -bj A ) v-y A 



(3) 



where P is the normalized line profile, M is a profile mask (M=l 
unless cosmic ray removal is used, see below), T is a fractional 
profile (0 <r< 1, T= 1 unless a cross-talk correction is applied, 
see below), and V' ^ = V j A . In this case the flux is integrated 
across j using a wider aperture than with the standard extraction 
(x w , cf. Sect. [2j. The variance of the modified optimal spectrum 
estimate is, moreover, 



Ka = 



(4) 



After fa is calculated using V' j A = V j A (Eq.[T) we instead use 

V'ojA = \fuPij*\ /gjA + r ) A (5) 

and iterate the solution «2 times. Again, following H86 the de- 
fault is «2 = 1. 

Next we describe how we remove cosmic rays and correct 
the integrated flux for cross-talk. Both methods require an in- 
spection of the resulting outcome, in order to see that it is sat- 
isfactory. These are therefore options that must be switched on 
separately in p3d. 

Removal of cosmic ray hits In order to remove cosmic ray hits 
we follow the approach of H86 and first iterate the integrated 
flux (Eq.O at most «cr times. In each iteration we mask at most 
one pixel (j) with the highest value which satisfies 



M m l(d jA - b jA - f u f - , 



V' „ > 



CR " oJA 



by setting M(j A = 0. Here <tqr is a threshold that defines how 
large the deviation of one single pixel must be to be classified as 
a cosmic ray hit. In comparison to H86 we found that it is nec- 
essary to use values larger than <tcr = 5. This method sometimes 
has difficulties removing cosmic rays if emission lines are very 
strong. Unless <xcr is set high enough (k. 10) pixels in the line 
center of profiles with high intensity may be removed as cosmic 
rays; resulting in a strong decrease of the integrated flux, as the 
lost flux is not replaced with any interpolated (or expected) flux. 
The default value on the number of iterations is «cr = 2. 

Correcting extracted flux when there is overlap between 
nearby spectra on the CCD In order to correct for cross-talk 
we iterate the spectrum extraction. In every iteration we first cal- 
culate a total profile across all m\ spectra of group /. Thereafter 
we calculate a fractional profile r,-^ that for every contributing 

pixel indicates which fraction, of the flux (dj A - bj,ij and the 
variance V' j A , belongs to line profile i. The spectrum extrac- 
tion is iterated at most «ct times, or until the maximum relative 
change of the calculated extracted flux of all spectra in group /, 



of two consecutive iterations, is < cr cx . The default number of it- 
erations and value of the threshold are «ct = 15 and <x CT = 10~ 5 . 
The number of required iterations depends on the data and the in- 
strument, but it appears that less than five iterations are required, 
typically, using data of the PPAK IFU. 

3.2.3. Multi-profile deconvolution optimal extraction (MPD) 

In addition to the modified optimal extraction p3d can also use 
the multi-profile deconvolution method of S10. In comparison 
to our method above all line intensities at one wavelength are 
here solved for simultaneously. Specifically, for every group of 
spectra the method is to minimize the residual 



Rf,A - 2 Z 



1 ri (<1.U-I> ;i -X, /ii/'//l) 2 



Voja 



Assuming dRf tA /dfj A =0 we then find and solve 

/ \ fiACijlA = bu 
i 

for the intensities fa, where 

p itfP m 



(6) 



(7) 



CfjlA 



-z 

m 



v, 



o.jA 



and bri\- > 

UA ^ y 



X0 



v ,jA 



Considering the variance we calculate it using a similar ap- 
proach by minimizing the residual 



Rv.a = 2 Z V ja - r~ jA - ^ V sJA P ijA 

/(!) V i 

Assuming dRyj/dV s j A =0 we then find and solve 



for the variances V s j A , where 

cvjia = 2 PijxPijA and b VJA = ^ (v o j A - rf jA } P 
jd) AO 



(8) 



(9) 



The other two spectrum extraction methods use a relatively 
limited size of the aperture. In this method we use the aperture 
width 

i w = 2n mp d<i/bin-|- + 1 , 

where « mp d is the number of neighbor profiles, on either side of 
every spectrum, which are considered. By default n mp d = 1, but 
it can be increased to include more spectra if line profiles are 
broad. The profiles P are calculated according to the description 
in Sect. [3X2] 

Following the approach of S10 the systems of equations 
(Eqs. [7] and O are solved for fa using a tri-diagonal solver if 
«mpd = 1- For n mp d > 1 there is a choice of solving the equa- 
tions using either a solver for a sparse diagonal matrix or singu- 
lar value decomposition. 
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3.3. Preparing a dispersion solution 

Depending on the optical path through instrument and fibers re- 
sulting spectra of separate spatial elements of an IFU are shifted 
and stretched relative to each other. An image with extracted 
spectra of an arc lamp exposure shows emission lines as curves 
across the cross-dispersion axis. In order to find a dispersion so- 
lution, which is unique to an IFS exposure, all spectra should be 
aligned and stretched to use the same starting wavelength and 
size of wavelength bins. 

In order to calculate a dispersion mask p3d requires a list 
of lines with known wavelengths, some information about the 
expected wavelength range of the used setup, and an arc lamp 
image. In brief, the list of lines is first modified to match pixel 
positions of corresponding entries in the data. Thereafter a poly- 
nomial is fitted between pixel positions and corresponding wave- 
lengths. Next we describe the individual steps in more detail. 

In a first step a list of lines, with well known wavelengths, is 
selected; either automatically using information in the arc image 
file header, or as defined by the user. For PMAS, SPIRAL, and 
other instruments where the canonical reflective grating equa- 
tion applies, an initial estimate of the wavelength range is cal- 
culated using information about the instrumentation setup in the 
data header. For VIRUS-P the wavelength range is pre-defined 
as 3620-5910AQ 

In a second step a line mask curvature is determined by cal- 
culating positions for one or two lines in the arc frame. The algo- 
rithm to do this involves locating the maximum position of one 
line in one spectrum, and then track this maximum through all 
other spectra. Emission lines of separate fibers in PMAS data us- 
ing the 2kx4k-CCD are, to first order, only shifted in wavelength 
relative to each other. In this case we found that it suffices to se- 
lect one emission line in the data in order to calculate the cur- 
vature accurately enough. VIRUS-P, SPIRAL, and PMAS using 
the 4kx4k-CCD are different, for these instruments the disper- 
sion varies across the IFU surface. In this case we selected two 
clearly separated emission lines in the data, in order to calculate 
both the curvature and the change of dispersion. 

We present examples of curved arc lines, in the blue wave- 
length range, for LARR in Fig. [4^ (for lines of mercury) and for 
VIRUS-P (bundle 2) in Fig. 2J; (for lines of cadmium and mer- 
cury). In order to illustrate the change of dispersion of VIRUS-P 
we, in this case, calculated the curvature using one line, instead 
of two. As Fig. |4j; shows the curvatures of lines in the line list 
badly matches the data in the redder part of the image. If two 
lines are used instead to calculate the curvature the match is ex- 
cellent across all spectra. 

In a third step the line mask is shifted along the dispersion 
axis in order to achieve a rough match between entries of the 
line mask and lines in the arc image. If lines cannot be matched 
across the entire dispersion axis the constant pre-estimated dis- 
persion can be changed manually. With the PMAS 4kx4k-CCD, 
in particular, the dispersion changes between the blue and the red 
ends of the detector. In this case it is necessary to fit a prelim- 
inary non-linear dispersion solution to an interactively selected 
set of lines. If the selected line mask, moreover, contains satu- 
rated lines, or if there are more entries in the input line mask than 
are visible in the data, such entries can be removed interactively 
in a fourth step. 

More precise pixel positions of entries in the line mask are 
then calculated by correlating every line in the line mask with the 



7 This wavelength range can, if necessary, easily be changed by pro- 
viding p3d with alternative lower and upper values. 



data, this is done for every spectrum. In this fifth step more pre- 
cise pixel positions are by default calculated using a Gaussian 
function fitting. Alternatively they can be determined using a 
much faster iterated average-weighting scheme. The region that 
is searched for an intensity maximum normally spans a range of 
i w pixel on the dispersion axis. If emission lines are too tightly 
packed, in either the line mask or in the data, this step may fail; 
in this case it is recommended to remove a few entries in the line 
mask. 

The sixth, and final, step is where a dispersion mask is cre- 
ated. A (linear) polynomial of pre-defined order p is fitted to the 
pixel positions at all wavelengths, for every spectrum, in the line 
mask. The polynomial order can be set to any value, the default 
is to use a low order with p = 4. Residuals of the fits are stored to 
the data reduction log file, and can also be viewed interactively. 
The dispersion mask is saved as p + 1 fitting parameters of every 
spectrum. 

Finally, we show an example of how an extracted arc lamp 
frame appears before and after the dispersion solution is applied 
in Figs. [4^ and Hfc. Note the extremely noisy line in the right 
part of Fig.|4}5. This line lies outside the range of selected lines, 
which were used in the creation of the dispersion mask, and was 
therefore inadequately calibrated. 

3.4. Flat fielding the extracted IFU data 

The optical path and transmission efficiency of individual spa- 
tial elements across an IFU typically vary. Differences appear 
both as wavelength-dependent variations and as a variance in 
the fiber-to-fiber throughput. For multiple-detector instruments, 
such as VIMOS, there is likely also a difference in the detector- 
to-detector throughput. After spectra are extracted and wave- 
length calibrated p3d can correct for these variations by normal- 
izing the data with an extracted flat field image, which combines 
the required corrections. 

A correction for a variable fiber-to-fiber throughput is ap- 
plied by dividing every spectrum of a flat field image with its 
mean spectrum. Correcting for wavelength-dependent variations 
every spectrum is at first smoothed across the dispersion axis us- 
ing a box-car of width wg pixel. Thereafter it is replaced with a 
(linear) polynomial fit of order pg. Finally, the flat field is nor- 
malized with the mean value of all elements. Only non-zero ele- 
ments are used with these operations. This smoothing minimizes 
the amount of noise that is added to the flat fielded spectra (see 
e.g. lBeckerll2002t iRoth et aTll2005l) . This smoothing can, if re- 
quired, be switched off (by setting ws = and pff = 0). 

Before a normalized flat field is created p3d, by default, first 
calculates a trace mask for this task using the same image. If twi- 
light exposures of flexure-affected instruments are used as flat 
fields this is an important aspect as it assures that the proper 
traces are used. The extracted flat field can, if required, be wave- 
length calibrated using a separate dispersion mask. 



3.5. Object extraction 

The only prerequisite to extract spectra of object data is a trace 
mask. Cosmic ray hits can, currently, be removed either by al- 
lowing p3d to combine a set of raw data images, or by using 
the cosmic ray removal option of the modified optimal extrac- 
tion procedure for separate images. Alternatively, for separate 
images they could be remov ed in advance , outsi de p3d, using, 
for example, the approa ch of Ivan Dokkuml J200ll who calls his 
routine LA. Cosmic) or lPvchl b200W . In principle a cosmic -ray 
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Fig. 4. This figure shows three extracted arc frames of: a) LARR, b) LARR after the dispersion solution was applied, and c) VIRUS- 
R Colors are modified to enhance various features. In each panel the vertical line with a disk indicates a preset reference spectrum 
and wavelength bin. In panel a the arrow points at the merging point of a shorter (Is; left part) and a longer (10s; right part) exposure, 
note that this straight line is curved in panel b (left arrow). In panel b the right arrow points at an emission line that was outside the 
calibration range. In panel c the left arrow indicates the curved shape of a line in the line mask, after the line geometry has been 
calculated. The two right arrows indicate a mismatch between the line mask and the data for redder wavelengths. For further details 
see Sect. [331 



removal option can also be added to the multi-profile deconvo- 
lution extraction procedure. We do not recommend to use this 
approach, however, since it is difficult to remove such hits effi- 
ciently and also keep all flux. Instead, for a future version of p3d 
we propose to use a mask (calculated using e.g. LA. Cosmic) 
that indicates cosmic -ray affected pixels. The intensity error of 
the masked pixels are set to a high value before the extraction, 
and will thereby be given a minimal weight using either of the 
optimal extraction routines. Errors are, furthermore, calculated 
and stored in a separate file for all intermediate products, but 
only if a master bias is specified. 

After their extraction spectra are wavelength calibrated if a 
dispersion mask is specified. If the object data contains sky emis- 
sion lines p3d optionally calculates a shift of the dispersion mask 
based on the known wavelengths of those lines. The offset is cal- 
culated by first fitting all present sky emission lines in all spectra 
with a Gaussian profile. The offset is then taken as the median 
of the difference between the center positions and the expected 
wavelengths. The error of wavelength calibrated spectra should, 
moreover, be considered a lower limit since pixel values are not 
cross-correlated when interpolating the dispersion solution to a 
common base wavelength. 

p3d stores all data in row-stacked-spectra (RSS) format, 
although final spect ra can optionally also be saved in the 
E3D-format (see e.g. iKissler-Patig et al.ll2004b . When the RSS- 
formatted file is used to view or analyze the data further, outside 
p3d, it is necessary to use a separate table specifying the posi- 
tions of the spatial elements. The tables of the p3D-distribution 
in this case provide all necessary information. If the data need to 
be corrected for effects of differential atmospheric refraction we 
recommended to make t his correction befor e the da ta is flux cal- 
ibrate d. The approach of lFilippenkol d!982l cf. e.g. ISandin et al.l 
120081) provides the most straightforward approach as it only re- 
quires information about observing conditions (the required in- 
formation is mostly found in the data header). 



4. Program validation 

In this section we present the outcome of our tests of p3d. 
Primarily we used simulated data since the outcome then is 
known. In our study of properties of the trace mask and the 
dispersion mask we also compared outcome of p3d with cor- 
responding outcome of iraf for data of the LARR IFU. A de- 
scription of the observational s etup f or used data, in this case, 
can be found in iRelaflo et al.l d2009t) . Th e iraf data-reduction 
follows the scheme that is presented by lAlonso-Herrero et al.l 




dispersion axis [px] dispersion axis [px] 

Fig. 6. One trace is shown for the LARR IFU in panel a, and for 
the PPAK IFU in panel b. The position on the cross-dispersion 
axis is in both cases drawn as a function of the position on the 
dispersion axis. The solid line shows the smoothed trace and the 
dash-triple-dotted line the non-smoothed trace. Compare with 
Figs.|5li&|5p, cf. Sect. 



d2009h . with some minor modifications (see IRelaflo et alj|2009l 
for details). 

At first we discuss the accuracy of the trace mask in Sect. 14. II 
and study the importance of using a cross-talk correction in 
Sect. 14.21 Thereafter we compare the outcome of the wavelength 
calibration of p3d with outcome which is created using iraf in- 
stead in Sect. 14.31 In Sect. 14.41 we evaluate the accuracy of the 
three spectrum extraction methods. Finally, in Sect. l4.5l we com- 
pare resulting spectra using the different spectrum extraction 
with observations of planetary nebulae. 

4.1. Estimating the accuracy of calculated trace masks 

Accurately determined traces, i.e. spectrum center positions, are 
necessary in order to extract spectra properly (cf. Sect. l4.4l ). The 
automatic tracing algorithm of p3d is mostly able to locate all 
spectra at one wavelength, without any interaction, in the first 
step. During the spectrum tracing across all wavelengths, of the 
second step, the accuracy of the result depends on how well the 
profile center positions are determined. 
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Fig. 5. This figure illustrates the influence of the cross-spectrum smoothing of step 2 of the tracing algorithm. Four trace residuals, 
which are calculated from typical continuum images, are shown for each IFU as a function of the dispersion axis: a) LARR, b) 
PPAK, c) SPIRAL and d) VIRUS-P; the traces are selected uniformly across the IFU surface. The black (light gray) solid lines 
show the respective trace where the iu a = 1 1 -smoothed (io a = 21-smoothed) trace mask is subtracted from the non-smoothed mask 
(w a = 1). In order to plot all four traces using the same axis they are offset by half a pixel from each other. Dark gray horizontal lines 
are guides. For further details see Sect. 14.11 



In Fig.|5]we compare calculated traces, which are smoothed 
across several spectra (w. d = 11, 21, cf. Sect. I3.l!2l , with non- 
smoothed traces (io a =1). The difference between the two sets 
are $ 0.2 pixel in every case, except at the blue and red ends 
where differences are up to about 0.5 pixel. All four IFUs show 
semi-periodic offsets in the non-smoothed traces, that is due to 
inaccurate pixel subsampling. The oscillation period of a trace 
corresponds to the number of illuminated pixels of a row on the 
detector. We illustrate this for one typical trace of the LARR IFU 
in Fig. and for one similarly typical trace of the PPAK IFU 
in Fig. [(J. The trace of the LARR IFU extends across about 10.5 
pixels on the cross-dispersion axis, corresponding to 10.5 oscil- 
lation periods, which are seen in the second trace from the top in 
Fig.|5^. Likewise the PPAK trace extends across about 2 pixels, 
that equals the number of oscillation periods of the two middle 
traces in Fig. |5j). In addition to offsets of inaccurate pixel sub- 
sampling an additional offset is caused by the smoothing across 
several spectra. In this context the exact value of w a seems unim- 
portant as differences are both lower and higher when the traces 
of either smoothed mask are compared with the non-smoothed 
mask. 

Next we compare traces created with p3d and iraf for the 
LARR IFU. At first we show four residual traces where the 
trace mask of p3d was subtracted from the trace mask of iraf, 
cf. Fig. [7^. As the figure shows these residuals are always 
< 0.2 pixel, and mostly < 0.1 pixel. The residuals of the non- 
smoothed trace mask show slightly higher variations. We also 
show residuals using fitted traces (dotted line; where the line 
center positions were fitted together with the intensities when 
calculating the profiles). In this case the fitted traces differ 



slightly from both the fixed traces and those of iraf. Note that the 
fitted traces do not oscillate across the dispersion axis. In com- 
parison to the outcome of the Gaussian weighting (solid lines 
in Fig. [7j this indicates a negligible influence on Gaussian-fitted 
center positions due to pixel subsampling. 

For this test we also fitted the traces of p3d with a fourth 
order linear polynomial, which we then subtracted from the trace 
mask of iraf. We show the result in Fig. |7j5. In this case the 
non-smoothed traces show a better agreement with the traces of 
iraf. As both panels in the figure show the difference is minute 
when the fitted traces are replaced with a polynomial fit. The 
minimum and maximum values of the residual of all spectra are 
-0.17 and 0.22 pixel. The mean and standard deviation are -0.043 
and 0.044 pixel. In Fig. [9] we show an image of the residuals of 
all spectra. The discrepancy is, again, the largest at the blue and 
the red ends. 

In order to compare the calculated center positions with the 
raw data we indicate the center positions of each approach in 
Fig-El for two wavelengths. We selected the bluemostpart of the 
second trace from the top in Fig.|7](lefthand side panel), and the 
redmost part of the topmost trace (righthand side panel), since 
this is where the difference is the largest in the shown residuals. 
In the former case the fitted center position of p3d lies closer 
to the position of iraf, the difference compared with the pre- 
calculated trace mask position is 0. 1 1 pixel. In the latter case the 
same difference is only 0.012pixel, while the position of iraf is 
more offset. The profile center positions of p3d are here more 
accurate than those of iraf. 

We conclude that the error of (fixed) traces, which are cal- 
culated using p3d, should be 5= 0.2 pixel. For spectra, which are 
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Fig. 7. Four trace residuals are shown for a typical continuum image of the LARR IFU as a function of the dispersion axis: a) 
trace(iRAF)-trace(p3D), b) trace(iRAF)-fit(trace(p3D)); the traces are selected uniformly across the IFU surface. The black (gray) lines 
show the residuals using the smoothed (non-smoothed) traces of p3d. The dotted line shows residuals of traces that were fitted anew 
during the profile calculation. In order to plot all four traces using the same axis they are offset by half a pixel from each other. The 
dark gray horizontal lines are guides. For further details see Sect. 14. II 
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Fig. 8. Profile center positions are shown across the cross-dispersion axis for iraf, p3d, and fitted profiles of p3d for two wavelength 
bins of two different spectra on the detector. The raw data is drawn as a histogram, and the Gaussian fitted profile is shown as a 
continuous curve. For further details see Sect. 14. II 
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Fig. 9. In this image we show the difference between a trace mask that is calculated using p3d and a corresponding trace mask that 
is calculated using iraf. For further details see Sect. I4.ll 



located close to the center of the detector, this value is likely bet- 
ter. At the blue end of the detector, for A 5= 4000 A the accuracy 
is likely worse due to a lower detector sensitivity. If high accu- 
racy is required, as may be the case when cross-talk is present 
(cf. Sect. I4.4b . the preferred method is to replace the fixed center 
positions of the trace mask with fitted positions. Note, however, 
that it may be difficult to achieve the highest accuracy of the 



spectrum extraction with flexure-affected instruments. Because, 
spectra might have moved on the detector during the time be- 
tween the calibration exposure and the object exposure. A solu- 
tion to this problem could be to re-center the calculated profiles 
using the object data. 
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Fig. 10. This figure illustrates effects of fiber-to-fiber cross-talk between neighbor spectra for typical idealized profiles of four IFUs: 
a) LARR, b) PPAK, c) SPIRAL (Blue arm), and d) VIRUS-P (bundle 2; bin-| =2). The abscissa shows a part of the cross-dispersion 
axis, while the ordinate shows the intensity in an arbitrary unit. The reference profile is drawn with a black solid line, and neighbor 
profiles with gray lines. In every case profiles, which are drawn with solid lines, have equal intensity. Profiles that are drawn with 
a dash-dotted (dash-dot-dot-dotted) line are 10 (100) times as intense. The summed spectrum is drawn with a thick light gray line. 
The center position and width of the aperture are indicated at the top of each panel. For further details compare with the values in 
Tableland also see Sect. 14.21 



Table 3. Parameters and outcome of our idealized study of the influence of cross-talk 



IFU 


bin t 




G 


nean 


G,- 


ngc 


-<v/bin t 


Si 


<5l/R 


<5l0/R 


<5l00/R 






[px] 


[px] 


[px] 


[px] 


[%] 


[%] 


[%] 


[%] 


LARR/V600/blue 


2 


2.00 


1.78 


±0.06 


1.6- 


-2.1 


7.0 


0.0004 


0.03 


0.13 


1.1 


PPAK/V600 


2 


2.25 


2.27 


±0.09 


2.1- 


-2.5 


5.5 


0.44 


2.6 


11 


110 


SPIRAL/Blue Arm a 


1 


2.30 


2.25 


±0.08 


2.1- 


-2.5 


4.0 


3.7 


3.8 


19 


170 


SPIRAL/Red Arm 


1 


2.30 


2.43 


±0.15 


2.1- 


-2.6 


4.0 


5.3 


5.6 


28 


260 


VIRUS-P/bundle 1 


1 


5.00 


4.73 


±0.47 


3.8- 


-5.4 


7.8 


5.2 


4.4 


23 


200 


VIRUS-P/bundle 2 


1 


4.20 


4.01 


±0.23 


3.7^.5 


7.6 


2.6 


1.4 


7.2 


67 


VIRUS-P/bundle 2 


2 


2.10 


2.07 


±0.10 


2.0- 


-2.3 


4.3 


1.5 


3.6 


18 


160 



11 Due to low intensities in the bluemost part we only used the redmost 1650 pixels on the dispersion axis. 



4.2. On the importance of correcting for fiber-to-fiber 
cross-talk 

The same calibration data that is used to calculate the trace mask 
is also used to calculate cross-dispersion line profiles, which are 
used in the optimal extraction. For most instruments there is 
some overlap between spectra on the cross-dispersion axis. The 
magnitude of the overlap depends on three factors: the spectrum 
separation (d), the width of the profiles, and the intensity of ev- 
ery spectrum. In order to estimate the influence of cross-talk we 
first calculated line profiles for different setups of each supported 
IFU. In this idealized study we assumed a Gaussian shape of the 
profile in every case. Note that this simplified treatment neglects 
the effect of extended profile wings, which are due to scattered 
light; typicall y at < 1% peak intensity, but large FWHM (cf. e.g. 
lBeckerll200l . 

Using data from different observing runs we first measured 
both the variation and the average value of the spectrum width 
across the IFU. We present the outcome in Table [3] The initial 
FWHM of the profiles of every group (Gi n j = G w /bin-| -) was taken 
from Table[2] The calculated average width and its standard devi- 
ation are given in Col. 4 (G me anX an d the full range of measured 
widths in Col. 5 (G range ). Note that the width varies across the 
detector for all IFSs, in most cases across both axes. This vari- 
ation should be kept in mind when interpreting the percentages 
we present next. 

For each IFU configuration we then calculated a set of 3-4 
line profiles using the spectrum width G mean and the spectrum 



separation d. We integrated the flux for one (reference) profile 
across an aperture of pre-defined width x'^ {x' w = x w /bin^, see 
Table |2]i and calculated the fraction of the flux which fell out- 
side the aperture (6{). Additionally, we calculated the fraction of 
increased flux inside the aperture, due to the two neighbor pro- 
files, in order to estimate the cross-talk contribution. Doing this 
we assumed that the intensity of the right-hand-side profile is 1, 
10 and 100 times higher (5i/r, 5io/r, 5ioo/r) than the reference 
profile. The outcome for each IFU setup is given in Table [3] We 
additionally illustrate the profiles of four setups in Fig. [TO] 

This study reveals several important results. For the LARR 
IFU it is evident that a cross-talk correction is unnecessary. Even 
if there is a strong intensity gradient across the IFU, with an in- 
tensity ratio of 100 between two neighbor spectra, the amount 
of cross-talk is only about 1%. The amount of flux outside the 
aperture of the reference spectrum of the PPAK/IFU is negli- 
gible (0.44%). Since fibers are more tightly packed, however, 
the amount of cross-talk is significant. If intensity gradients are 
moderate (and ratios are < 10) the intensity of the reference pro- 
file is increased by about 12% due to cross-talk. With the remain- 
ing IFU configurations the amount of flux outside the aperture 
is 1-5%, while the cross-talk contribution to the reference pro- 
file increases to 7-28% for moderate overlap, and to 67-260% 
for strongly overlapping profiles (intensity ratio of 100). With 
the exception of the LARR IFU it is clear that it is necessary to 
correct for cross-talk in order to achieve any level of accuracy in 
integrated fluxes of weaker regions of the object. Although, even 
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for the LARR IFU an optimal extraction could be advantageous 
if it is important to extract weak lines accurately (cf. Sect. 14.5b . 

I Arlington- Smith & Content! (1 19981) argue that a cross-talk 
correction is unneccessary when neighboring spectra on the 
detector correspond to well-sampled neighboring positions on 
the sky, as is the case with the LARR and the SPIRAL 
IFUs (and also the VIMOS and the FLAMES -ARGUS IFUs). 
This argument holds particularly well whenever precision spa- 
tial sampling is not the driver for a given application; us- 
ing, for example, adaptive binning with SAURON observations 
dCappellari & Copinll2003l) . However, owing to properties of a 
datacube, which can be seen as a stack of hundreds of images, 
IFS with proper spatial sampling has the potential of deliver- 
ing extremely accurate astrometry with milli-arcsec centroid- 
ing precision for point sources from ground-based observations 
dRoth et al.l2004l) . For such applications, where an accurate defi- 
nition of the point-spread-function is essential, cross-talk correc- 
tions are very important (ibid. Fig. 4). For the SPIRAL IFU S10 
(see Fig. 7 therein) also find unwanted artefacts in their extracted 
spectra when they do not apply a cross-talk correction. 

While spectra of most IFUs are arranged next to each other 
on both the sky and on the detector, this is not the case for the 
PPAK IFU. Nearby spectra on the sky are on PPAK placed at 
different locations on the CCD. When seeing prevents extreme 
intensity gradients for IFUs with densely packed elements, in- 
tensity ratios of as much as 1000 are easily achieved with PPAK 
data; observing, for example, the center of a planetary nebula 
and its thousand times fainter halo. In this case it is necessary 
to use a two-dimensional approach of data reduction in order 
to ext ract the weak component properly (cf. iBolton & Schlegej 
120091) . 

4.3. Comparing the wavelength calibration correction 
between p3d and iraf 

Here we compare how well the wavelength calibration correc- 
tion of p3d matches a correction that is carried out with iraf 
instead. Using data of the LARR IFU we calculated three dis- 
persion masks, which we then applied to the extracted spectra of 
the respective arc image. We used the same line list with both 
tools. With p3d we used a fifth order linear polynomial, and with 
iraf a fifth order Legendre polynomial. The observations were 
done using the V300 grating of PMAS, where the dispersion is 
1.67 A/pixel. 

We evaluated the accuracy of either correction by fitting 
a Gaussian profile to all spectra, using several arc lines in 
the extracted and wavelength-calibrated images. Specifically 
we fitted two mercury lines with the rest wavelengths 
A = 4358.328, 5460.735 A, and two neon lines with A = 
5975.534, 6598.953 A. The standard deviation of the center po- 
sitions were found to be similar for the two tools; they vary 
between 0.06-0.31 A (iraf) and 0.09-0. 13 A (p3d) for the four 
lines and all spectra of the three images. In Fig.QT]we plot his- 
tograms of the fitted wavelengths for one of the dispersion masks 
and Hg /15461 and Ne /15976. For the mercury line (Figs. [TO and 
fTTb ) the respective tool provides the average (and standard devi- 
ation) value of 5460.71 (0.15, iraf) and 5460.69 (0.1 1, p3d). The 
results are evidently very similar in this case. For the neon line 
the corresponding values are 5975.43 (0.07, iraf) and 5975.60 
(0.12, p3d). In this case p3d shows a larger scatter of the values. 
The residual plots (Figs. [TO: andQjJ) confirm that the standard 
deviation of values of iraf are systematically smaller than with 
p3d, but differences are small. The wavelength-calibration cor- 



rection of p3d can likely be further improved if additional care is 
taken to enhance the resampling algorithm to a common wave- 
length. 

The accuracy that can be achieved with the wavelength cal- 
ibration depends on several factors. Four such factors are: the 
spectral resolution, the number of entries in the line mask and the 
accuracy of the center pixel positions, the fitting function (and its 
order if it is a polynomial), and properties of the final interpola- 
tion to a common wavelength for all spectra. Selecting a line list 
we have the following four recommendations: lines must be iso- 
lated, lines with high S/N are preferred, lines that are (close to) 
saturated should be avoided, and lines must be distributed across 
the full spectral range of interest. 

4.4. About the accuracy of spectra which are extracted using 
the three different methods 

In order to measure the accuracy of the spectrum extraction we 
test the three extraction methods with idealized simulated data. 
With this motive we assumed that the line profiles are Gaussian. 
All spectra are perfectly aligned with the dispersion axis, and the 
intensity is invariant with wavelength. Hereby we only model 
one wavelength bin. 

Setting up our simulation we first defined a set of ten 
Gaussian profiles for each of the four IFUs. We used the av- 
erage instrument-specific profile widths, G mea n, which we cal- 
culated in Sect. 14.21 (see Table 0). The Gaussian profiles are, 
moreover, separated by d/bin-j pixels (see Table [2]). We set the 
intensities of all, but two, profiles to the same value; the third 
and sixth profile intensities are set ten and hundred times higher. 
Thereafter the profiles are scaled to a pre-defined value of the 
signal-to-noise (S/N) of the weaker profiles, and are summed up 
to create one spectrum. We refer to the scaled intensities as mod- 
eled intensities below. For the noise model we used a Poissonian 
noise distribution and the instrument-specific readout noise. In 
all subsequent measurements we used the intensity of the sec- 
ond, fifth, and ninth profiles. These profiles correspond to mea- 
suring a weaker line that lies next to a line that is 10, 100, and 1 
times as strong. 

Evaluating the simulation we used an approach that is similar 
to that of S 10. For each IFU setup we fitted three different sets of 
Gaussian profiles to the simulated data. For the first set we used 
the already known center positions and profile widths, and scaled 
the input intensities to achieve a S/N of 2-300. For the second set 
we used the known profile width and set the input intensities to 
correspond to S/N=100, we also introduced an error to the pre- 
determined profile positions of 0.0-0.7 pixel. In the third set we 
used the known profile center position and again used S/N=100. 
This time we introduced an error to the profile width, which are 
used in the fitting, of -0.4-0.6 pixel; a negative error corresponds 
to a narrower profile and a positive error to a wider profile. For 
every configuration we fitted the profiles 100 times using the 
p3D-routine of the respective extraction method, using as many 
realizations of the noise model, and saved the average intensities. 

In the following text we refer to the modified optimal ex- 
traction method as MOX (Sect. [3~X2b . We, likewise, refer to the 
multi-profile deconvolution optimal extraction method as MPD 
(Sect. [TO]). 

We begin our analysis with the PPAK IFU. In Fig. [T2h we 
show the result of the high-intensity simulations where we in- 
troduced an error to the pre-determined center positions of the 
three profiles. Using the aperture extraction it is seen that the ex- 
tracted flux, due to cross-talk, always is higher than the model 
input intensity. Note that the values at error 0.0 agree well with 
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Fig. 11. This figure shows histograms and residuals of fitted emission lines from two extracted and wavelength-calibrated arc images 
of LARR data. In the top (bottom) row we show plots of Hg /I5461 (Ne /15976). Panels a and d (b and e) show histograms of fitted 
wavelengths using p3d (iraf), and panels c and f show residuals between the fitted wavelength and A^. The residuals are shown for 
p3d (+) and iraf (x, these residuals are offset by 0.5 A). For further details see Sect. 14.31 
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Fig. 12. This figure shows measured-to-modeled flux ratios for the PPAK IFU (using bin) =2). The three panels show the flux ratios 
as a function of: a) the assumed error of the profile center position, b) the assumed error of the profile FWHM and c) the signal-to- 
noise. We show flux ratios for three sets of profiles where the neighbor profile intensity is stronger by a factor: 1 (gray symbols), 
10 (black filled symbols), and 100 (open symbols, for clarity these values are drawn without connecting lines; the exception are 
the aperture extraction flux ratios which are drawn with a dotted line, these are also offset with -1.0 to keep the values in the same 
plot). Flux ratios, which are calculated using the standard aperture extraction, are indicated with circles, and ratios of MOX (MPD) 
with triangles (squares). All values in panels a and b are calculated using S/N=100. The horizontal lines at the flux ratio 1, and the 
vertical line at error 0.0 in panel b, are guides. For further details see Sect. 14.41 



the values in Sect. 14.21 (see Table [3j we we did not use discrete 
pixels in that study). The measured flux is ^ 100% higher than 
the model flux with the highest intensity ratio. If there are strong 
intensity gradients in the data, and it is important to measure 
weak regions accurately, it is not recommended to use the aper- 
ture extraction method. Using MOX the accuracy is higher than 



with the aperture extraction. Although with MOX the error in- 
creases rapidly with the error of the pre-determined center posi- 
tion for the two profiles, which lie next to the intenser profiles. 
MPD shows the best performance of all methods at small er- 
rors. Although, for larger errors (^ 0.2 px) it is outperformed 
by the other two methods when the intensity ratio is near unity. 
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With MPD it is enough to keep the profile center error smaller 
than 0.2 pixel in order to achieve an intensity error that is smaller 
than 5% for all intensity ratios. Using MOX the corresponding 
error of the pre-determined center position must, in this case, be 
smaller than about 0.03 pixel. 

We study the influence of an accurate pre-determined pro- 
file FWHM in Fig. [TZb . Both optimal extraction methods show 
a similar dependence for positive errors. MPD is seen to give 
a slightly higher accuracy at reasonably low positive errors 

0. 3 pixel). When the pre-determined FWHM is too narrow, and 
the error consequently is negative, MPD shows smaller errors in 
the flux than for positive FWHM errors. With MOX the resulting 
flux errors are similar as for positive FWHM errors. In order to 
calculate fluxes, which errors are smaller than 5%, it is necessary 
to keep the FWHM error smaller than 0.05 pixel (0.2 pixel) with 
MOX (MPD). 

The dependence of the resulting flux error on the S/N is 
shown in Fig. [TZb (note that the line of the factor- 100 inten- 
sity ratio is offset by -1.0 on the ordinate). It is seen that val- 
ues changes little with the S/N, with the exception of low values 
(S/N£ 4). Note that for S/N^ 30, and the highest intensity ratio 
(100), the resulting fluxes, which are calculated using MOX, are 
about 3% lower than corresponding values, which are instead 
calculated using MPD. 

In Sect. l4.2l we showed that there is no need to correct LARR 
data for cross-talk even for high inter-profile intensity ratios, as 
spectra are very well separated. In order to provide a complete 
study we, nevertheless, test the outcome when such a correction 
is made anyways. 

Figure [T3h shows that MOX gives accurate results for small 
center position errors and higher intensity gradients. Although 
when the intensity ratio is unity MOX gives unreliable fluxes, 
which are ;£ 98% of the true flux for small errors, see Fig. [T3~b . 
The offset dash-triple-dotted lines (at flux ratio 1.4) show the 
calculated fluxes using MOX, without a cross-talk correction - 

1. e. for regular optimal extraction. The lines illustrate that flux 
errors are small, even at relatively high errors of the center po- 
sitions and the FWHM. It does not make any sense to use MPD 
or MOX with cross-talk correction with LARR data. In this case 
noise in the profile wings increase the errors of the measured 
flux, instead of decreasing them. 

With the SPIRAL IFU and the VIRUS-P IFU, see Figs.QjJI- 
[T3l . we see a similar behavior as for the PPAK IFU. With 
these IFUs and MPD it is necessary to keep errors of the 
pre-determined center positions ;S 0.02 pixel (SPIRAL) and S5 
0.15 pixel (VIRUS-P) in order to delimit flux errors to 5%. With 
MOX the center error should be ^ 0.01 pixel (SPIRAL) and 
< 0.03 pixel (VIRUS-P). The required precision of center posi- 
tions is very high for SPIRAL, with both methods, which is why 
it seems unrealistic to achieve this high accuracy in fluxes. In 
order to delimit errors we recommend to fit the spectrum center 
positions, along with the intensities, when calculating line pro- 
files for this IFU (cf. Sect. 14. 11 1. Both methods, moreover, require 
an accuracy in the profile width of ;£ 0.05 pixel in order to con- 
strain the measured flux error to 5%. The S/N-test for VIRUS-P 
(Fig-US) shows that the error of MPD is about half that of MOX 
for S/N;S 10. For the SPIRAL IFU significant flux errors (<4%) 
are only introduced, as a function of S/N, with MOX (Fig.Qjf)- 

We conclude that MPD is the preferred method of spec- 
trum extraction whenever cross-talk is present. The cases where 
MOX was found to outperform MPD can likely be explained by 
the smaller extraction width that is used with MOX (compare 
3c w with l w that were both defined in Sect. 13.21 ). Regardless of 
the chosen method of optimal extraction it is always important 
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Fig. 15. We plot parts of the bias-subtracted raw data, which we 
used to extract the spectra of Fig. [14] The counts are plotted 
for one wavelength across the cross-dispersion axis for data of 
LARR, panel a, and PPAK, panel b. In panel a the group of 
profiles drawn in gray is shown for the center wavelength of Hy; 
the black line is drawn for the corresponding wavelength bin 
5 pixels towards the blue of the center of Hy. The dark gray 
horizontal line marks the zero-level. In panel b similar profiles 
are shown about the center wavelength of Ha. Positions of four 
apertures are shown with gray indicators. For further details see 
Sect. 1431 



to minimize errors of the pre-determined center positions and 
FWHM. Errors of measured fluxes grow fast if there are signifi- 
cant intensity gradients in the data. We found that MOX is highly 
sensitive to accurate spectrum center positions (i.e. traces). MOX 
also introduces some error to data with low S/N. Although with 
LARR data, which is free of cross-talk, MOX is the preferred 
method (because of the way the extraction width of both meth- 
ods is defined within p3d), without the optional cross-talk cor- 
rection. Finally, with the exception of LARR data, aperture ex- 
traction always introduces significant errors to calculated fluxes 
of data, which contains any spatial intensity gradients. 

4.5. Comparing resulting emission line spectra of planetary 
nebulae using the three different extraction methods 

Finally, we compare the outcome of spectra, which were 
extracted using the different methods. In Fig. [14] we show 
wavelength-calibrated and flat-fielded spectra of data of LARR 
and PPAK from the central regions of two planetary nebulae 
(NGC 2392 fo r PPAK, and M 2-2 for LARR, which outcome 
is presented in Sa ndin et al.l l2008i). In the first case we reduced 
the LARR data using MOX (without correcting for cross-talk) 
and aperture extraction. As Fig. [T4k shows the spectrum of a 
single spatial element is less noisy when MOX is used instead 
of aperture extraction. The error bars of MOX are half as high 
than with aperture extraction. In the spectrum that is averaged 
over 12 spatial elements we see that the continuum level of the 
MOX-spectrum is lower than that of the aperture-extracted spec- 
trum (by about 55%). In this emission line object scattered light 
becomes significant for the weak continuum that is dominated 
by readout noise (compare the gray and black lines in Fig.[T5ji). 
Using aperture extraction more of the scattered light is included 
in the resulting spectrum, and weak emission lines are underes- 
timated. Note that both methods calculate the same intensity for 
Hy. 
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Fig. 13. Figure[12]is repeated for LARR (top row, panels a-c; bin t = 2), SPIRAL (middle row, panels d-f; bin t = 1) and VIRUS- 
P/bundle 2 (bottom row, panels g-i; bint = 2). All values in the three panels of the leftmost column are shown as a function of the 
assumed error of the pre-determined profile center position. Likewise, values in the three panels of the middle (rightmost) column 
are shown as a function of the assumed error in the pre-determined FWHM (signal-to-noise). In panels a and b we additionally plot 
flux ratios for MOX (dash-triple-dotted lines which are offset by 0.4 pixel), for optimal extraction, without cross-talk correction. 
The dotted lines in panels d and f (g and i) are offset by -3.5 (-2.5). 
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Fig. 14. Comparison of spectra, which are calculated using the aperture extraction, MOX and MPD. In panel a we show spectra of 
LARR and in panel b spectra of PPAK. The flux is shown as a function of wavelength - note that the two panels are plotted for 
different wavelength ranges. Solid black (gray) lines show spectra using MOX (aperture extraction). In panel b we used MPD to 
calculate the spectrum that is drawn with a dotted line. In panel a the two lower (upper) spectra use the left (right) y-axis. Error bars 
are shown for parts of the spectra in both panels. Wavelengths of a few emission lines (cosmic ray hits) are indicated with a vertical 
line and the name of the line (CR). For further details see Sect. 14.51 



In order to illustrate the result of cross-talk in PPAK data 
we show Fig. [T4b . In this case MOX and MPD give very sim- 
ilar solutions, except at the location of the cosmic ray hit, at 
A =i 6474 A. The flux at the wavelength of the cosmic ray hit 
was removed by the MOX-algorithm. The shown spectrum lies 
in a weak part of the central nebula on the sky, but on the CCD 
the neighbor spectrum lies in a more intense part of the nebula; 
compare the apertures 9 and 8 in Fig. [T5b . Using aperture ex- 
traction flux of the neighbor spectrum falls within the aperture 
of the weak spectrum, and increases the flux of Ha significantly 
above its true value. In this case we attribute the difference in 
the continuum levels of the aperture and MOX/MPD methods to 
both scattered light and cross-talk. The continuum-region error 
bars of MOX are on average 0.88 times as high as the aperture 
extraction error bars. 

5. Conclusions 

We have presented a new and general data-reduction tool for 
fiber-fed IFUs. Our goal has been to write a user-friendly tool 
that works stably with different sets of instruments and data. In 
comparison to similar tools a strong advantage of p3d is that it 
can find and trace all spectra on the detector, mostly without any 
user interaction at all. Using the methods of optimal spectrum 
extraction data of all IFUs can also be corrected for cross-talk, 
which arises due to overlapping spectra on the detector. Since 
the same procedures are used with all implemented IFUs it is, 
moreover, a straightforward task to compare the outcome of dif- 
ferent observations. Although p3d is based on the proprietary 
Interactive Data Language (IDL) its use requires no IDL license. 



All components of p3d work with all platforms supported by 
IDL. The program code can be downloaded from the project web 
site at http://p3d.sourceforge.net. 

In order to validate the program code we have tested the dif- 
ferent parts using both simulated data and corresponding out- 
come of iraf. We found that p3d produces results comparable to 
iraf. For all IFUs, with the exception of the lens array of PMAS, 
p3d is able to extract more accurate values than iraf since p3d 
can correct for cross-talk. 

Although p3d has so far been configured for four IFUs it is a 
straightforward task to extend it to work with additional instru- 
ments. If the new IFU is similar to the already implemented ones 
it is just a matter of setting up another set of instrument-specific 
parameters. The same concerns the level of functionality. p3d 
can with relatively small effort be extended to also handle, for 
example, flux calibration, correction for differential atmospheric 
refraction (for those IFUs where it makes sense), removal of cos- 
mic rays in individual images, sky subtraction, and account for 
scattered light. Most parts of p3d are, moreover, fast and only 
require on the order of seconds to execute on a typical worksta- 
tion. The calculation of line profiles and the optimal extraction 
algorithm, however, are more computationally intensive and re- 
quire on the order of a few minutes. The code execution time 
could in this case be shortened by moving the relevant parts of 
the IDL-code to compiled (and dynamically loaded) C-code. If, 
and when, these and other improvements will be implemented 
depends on the need and the interest of the community. 
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